!***************************************
finish
/clear
/filn, E321  
/title, vortex
/units,si	!ùʵλ
! ********* 趨*********
R1=775		! Ͳڰ뾶
t1=100		! Ͳ
R2=800		! ڰ뾶
t2=48		! Ǻ
L=95		! Ͳ߳                 
Lc=1200	        ! Ͳ峤
p=16 		! ѹ
E=2e5		! ϵģ
nu=0.3		! ϲɱ
nt=5		! ȷʷ
ns=30		! Ǿʷ
nc=30		! Ͳʷ
nl=5		! ɶʷ
ra=0.6		! ʷֱ
! *********ǰ********* 
/prep7
et,1,82			! 嵥Ԫ 
keyopt,1,3,1		! 趨Գѡ 
mp,ex,1,E		!ϵĵģ
mp,nuxy,1,nu		!ϵĲɱ
cyl4,,,R2,0,R2+t2,90	!ǲ
BLC4,R1,0,t1,-Lc+L	!Ͳ岿
wprot,,-90		!תƽ
wpoff,,,L		!ƶƽ
asel,s,loc,y,0,R2+t2	!ѡǲ
asbw,all		!ùƽָǲ
asel,s,loc,y,0,L	!ѡӦɾǲ
adele,all,,,1		! ɾӦɾĲ
k1=kp(R1,0,0)           !ȡɶεĽڵ1
k2=kp(R1+t1,0,0)	!ȡɶεĽڵ2
ksel,s,loc,y,L	
*get,k3,kp,,num,min	!ȡɶεĽڵ3
*get,k4,kp,,num,max	!ȡɶεĽڵ4
Alls		        !ȫѡ
a,k1,k2,k3,k4		!ɹɶ
!  ʷ
lsel,s,loc,y,-Lc+L      !ѡںص߶
lsel,a,loc,y,0
lsel,a,loc,y,L
lsel,a,loc,x,0
lesize,all,,,nt,,,,,1   !趨ںʷ
lsel,s,length,,Lc-L	!ѡͲڱھ߶
lsel,r,loc,x,R1
lesize,all,,,nc,1/ra,,,,1!趨Ͳڱھ߶εʷ
lsel,s,length,,Lc-L	 !ѡͲھ߶
lsel,r,loc,x,R1+t1
lesize,all,,,nc,ra,,,,1		!趨Ͳھ߶εʷ
lsel,s,radius,,R2		!ѡھ߶
lsel,a,radius,,R2+t2
lesize,all,,,ns,ra,,,,1         !趨ھ߶ʷ
lsel,s,loc,y,0,L		!ѡɶھ߶
lsel,u,loc,y,0
lsel,u,loc,y,L
lesize,all,,,nl,,,,,1 		!趨ɶھ߶ʷ
alls
amesh,all		        !ʷ
fini
! ********** **************
/solu
nsel,s,loc,y,-Lc+L		!ѡͲ׶˸ڵ
d,all,uy		        !ʩλԼ
nsel,s,loc,x,0		        !ѡǶԳϵĸڵ
d,all,ux		        !ʩˮƽλԼ
lsel,s,line,,3,8,5
lsel,a,line,,10
nsll,s,1
sf,all,pres,p		        !ʩѹغ
alls
solve		                !
fini
!**********  **************
/post1
PLNSOL, S,EQV, 2,1.0		!ʾӦͼ
fini
